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We examine the dynamics of a correlated quantum dot in the mixed valence regime. We perform 
numerically exact calculations of the current after a quantum quench from equilibrium by rapidly 
applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibra¬ 
tion times and saturates upon the decrease of temperature at all times, indicating Kondo behavior 
both in the transient regime and in steady state. The time-dependent current saturation temper¬ 
ature matches the Kondo temperature at small times or small voltages; a substantially increased 
value is observed outside of linear response. These signatures are directly observable by experiments 
in the time-domain. 

PACS numbers: 73.63.Kv,72.15.Qm,02.70.Ss,05.60.Gg 


The Kondo effect is an intrinsically many-body phe¬ 
nomenon in which localized and itinerant electrons form 
a strongly correlated state m which shows signatures 
in thermodynamic, spectral, and linear response trans¬ 
port properties Originally introduced to describe 

the low-temperature properties of magnetic impurities 
embedded in a non-magnetic bulk [1], it has since been 
observed in a wide variety of correlated electron systems 
ranging from impurities adsorbed on surfaces [3 |7j to 
molecular transistors [5] , heavy fermion materials [5] , and 
mesoscopic quantum wires uni- Kondo behavior is also 
observed in semiconducting quantum dot heterostruc¬ 
tures, where a confined interacting region is coupled to 
non-interacting leads [IIHIl]- Gating of these systems 
allows to study Kondo physics over a wide parameter 
range [TT]. 

Kondo correlations only emerge at low temperature. 
Their onset is characterized by the Kondo temperature 
Tk that can be defined as the temperature at which the 
zero-bias conductance reaches half of its low-temperature 
value [Ij. In thermal equilibrium, the Kondo problem is 
well understood and quantitatively described by a range 
of analytical and numerical methods 

When a Kondo system is driven out of equilibrium, ad¬ 
ditional phenomena appear. For instance, the application 
of a (time-independent) bias voltage splits the Kondo 
peak [HIMI] and shows signatures in the double oc¬ 
cupancy and magnetization [32] ■ At large voltages the 
voltage dependence of the conductance decreases on an 
energy scale comparable to Tk [33H36] and its temper¬ 
ature dependence saturates at temperatures above 
[37H3S|. 

How observables evolve in time after a rapid change of 
parameters gnmi] and how they decay to their steady 
state limit is an open question. Recent experimental 
progress in the measurement of time-dependent quan¬ 
tities on ever faster time scales has enabled experimen¬ 
tal studies of such transient dynamics [42ll45] . making a 
theoretical description of quenches from correlated initial 
states important. 


In this paper we provide numerically exact results for 
the real time evolution of a Kondo system following a bias 
voltage quench from a correlated equilibrium ensemble to 
a non-equilibrium steady state. We show results for tran¬ 
sient and steady state currents and populations at tem¬ 
peratures ranging from T ^ Tk to T ^ Tk, which have 
not previously been accessible in numerical calculations. 
Our results, enabled by recent advances in numerically 
exact QMC methods, illustrate the time evolution of the 
current and its saturation temperature from the equilib¬ 
rium Kondo temperature to the increased steady state 
value and predict experimentally observable signatures 
of many-body correlations far from linear response. 

Quenches of initially uncorrelated systems have been 
examined by a variety of state of the art theoretical meth¬ 
ods. In such systems, electronic correlations are gradu¬ 
ally established as time progresses. Studies by the it¬ 
erative path integral approach (ISPI) [46ll5Ij . real-time 
renormalization group |52H54] . hierarchical equations of 
motion (HEOM) [55H37] , flow-equation methods |SS| , ex¬ 
act solutions in solvable limits (53] |3D| , perturbation the¬ 
ory m, time-dependent Gutzwiller approach [f>2L Ibd] . 
DMRG im [65], NRG |66|, continuous time quantum 
Monte Garlo (QMG) [5311371170] and bold-QMG [5711711 - 
m] have shown the dynamical build-up of Kondo cor¬ 
relations on exponentially long time scales, the time- 
dependent Kondo cloud formation in the leads [zano], 
and characterized in detail the steady state properties, in¬ 
cluding the current-voltage characteristics, voltage-split 
spectral functions and temperature dependence of the 
conductance. 

Quenches from a strongly correlated initial state pose 
a greater challenge, as the initial solution of an equi¬ 
librium many-body problem is required. At high tem¬ 
perature T > Tk, the dynamics following the voltage 
quench has been described by QMG [77[ [75] and HEOM 
[sa [73150]. These simulations have demonstrated time- 
dependent currents and spectral functions and related 
current oscillations to the applied voltage [5T]. At T = 0 
these results are complemented by DMRG [53H57] and 
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FIG. 1. (a) Interacting quantum dot embedded into a tun¬ 
nel junction. The dot is described by a level spacing Ed and 
a charging energy U. The leads are described by the half¬ 
bandwidth D and the chemical potential and Ed = 0. (b) 
At time t = 0 a voltage V is instantaneously applied to the 
system, changing the bath Fermi level to i:V/2. 

NRG [HSl IH3 IHS] simulations, recently extended to fi¬ 
nite temperature |87j . which show transient dynamics 
between the ground and steady states. Semi-analytic 
techniques, including NCA and OCA and the 

time-dependent Gutzwiller approach [62l [63] provide ap¬ 
proximate results for all values of voltage and tempera¬ 
ture. Their precision requires careful assessment jSH |S3| . 
In principle, QMC simulations can cover the full param¬ 
eter range of the system. However, so far these calcula¬ 
tions were limited to high temperature and short times 
by sign problems. We have overcome the first limitation 
in this work by successively normalizing to a sequence of 
reference systems at progressively higher expansion or¬ 
der, allowing us to gradually reach temperatures an order 
of magnitude below Tk- 

Model. We consider a quantum dot attached to two 
metallic leads, schematically shown in Fig. Both the 
quantum dot and the leads are initially in thermal equi¬ 
librium at temperature T. We describe the system by an 
Anderson impurity model 

H{t) -I- —V{t))cl^i^^Caka + Ht + Hjj (la) 

ak(T 

Ht + dicaka) (lb) 

k(7 

Hd ='^eddld„+ Un^ni. (Ic) 


the dot and leads is described by the effective coupling 
parameter F = 7rV^/(2ZI) and the hybridization function 

A^(t, t') = ^ Aa,j(t, t') = ^ f dojA{uj) 

a a ^ 

X [6»c(i,i') -/M]exp dt"{uj - . 

( 2 ) 

denotes the integral over the L-shaped Keldysh con¬ 
tour with an imaginary branch added to take into ac¬ 
count the initial correlated conditions [M], and /(w) = 
(exp(a;/r) -b 1)“^ is the Fermi function. 

We describe the evolution of the system by a hy¬ 
bridization expansion in the tunneling term Ht, which 
we stochastically sum to all orders using diagrammatic 
continuous-time Monte Carlo (CT-HYB) on the full 
Keldysh contour milSHlITH]. We perform an expan¬ 
sion for the partition function of the problem Z[IA\ = 
where every perturbation order n of 
expansion in Ht is characterized by conhgurations (fin 
consisting of 2n operators distributed on the contour and 
the impurity state I'F). Summing the series stochastically 
allows one to extract the time-dependence of the quan¬ 
tum dot charge [Zg as iV = Y. 4 ,{'^Ws^\'S))/{s), where 
J 24 . = is the Monte 

Carlo sign of every configuration and (s) = is the 

average sign. The sign (s) is equal to 1 in equilibrium, 
but decays exponentially with real time. In practice we 
resolve times of order of F”^. Calculations at larger times 
are exponentially more expensive. 

The main observable of interest is the current pass¬ 
ing through the dot. The QMC procedure should be 
modified for current calculations by replacing A with 
At(ti,t 2 ) = - St^d) + ^+hiA 2 )Sti,t and sam¬ 

pling the imaginary part of the weight of every configu¬ 
ration, yielding = Imw^^/jlmw^^ | [68|. From this we 
extract the quantity ZI{t) = 2lmZ[At] and, by dividing 
by Z, the current. 

In order to evaluate the current separate measurements 
of Z and ZI are required. These quantities are not di¬ 
rectly accessible by QMC, but are instead inferred from 
a normalization procedure, yielding 


c and d label electrons in the leads and on the dot, re¬ 
spectively, cr is the spin index, a = ±1 labels the left (-I-) 
and right (—) lead. V is the tunneling matrix element 
between impurity and leads, Sk is the lead dispersion, 
Ed the impurity level spacing, and U the impurity elec¬ 
tronic repulsion strength (charging energy). The chemi¬ 
cal potential is set to 0 at equilibrium, the instantaneous 
(quenched) application of voltage V{t) = V9{t) at time 
t = 0 changes it to a^. 

The leads are non-interacting and are described by a 
flat density of states A{uj) with half-bandwidth D and a 
smooth cutoff at the band edges [77]. The coupling of 


m 


r zi(t) 1 

■^ref 

, ■^ref-^ref (0 _ 





(3) 


where and /ref are the partition function and the 
current of a reference system, which are obtained sepa¬ 
rately. Each of the fractions in Eq. ([^ is obtained in a 
CT-HYB calculation as 

^ref _ E^R-eWref/lReW^I ^ Z ImWref / jlmw,^ | 

(4) 

The choice of the reference system is important: small 
ratios of /ref// and .Zref/.^ results in large relative error 
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bars. As temperature is decreased, the average pertur¬ 
bation order increases oc T“^, and the normalization to 
the first order of expansion in Ht |68j becomes unreli¬ 
able. We improve the algorithm in several ways: first, at 
high temperatures (perturbation orders < 10) the refer¬ 
ence system is a finite order NCA/OCA expansion. At 
lower temperatures (larger pertubation orders) we per¬ 
form a series of calculations, progressively increasing the 
perturbation order until convergence is achieved, yielding 


ZI = 


■ OO 

ZIm 

ZIm+{n-l)& 

Z/nca. 


ZInca, 


(5) 


with the order of the NCA expansion m = 7 and S as 
either 2 or 3 to maintain ratios > 10“^ at every step. We 
emphasize that our results are exact within the stochastic 
error bars for any choice of reference system; the choice 
affects only the computational cost. 

Results. We present results for the dynamics of a quan¬ 
tum dot in the mixed valence regime, Sd = 0. The Kondo 
temperature Tk is on the order of the coupling T El, 
and observables are expected to equilibrate within times 
oc accessible by QMC calculations [Ml[111172]. For 
convenient comparison to experiment we set the unit of 
energy close to 1 meV. Time is then given in meV“^ « 4 
ps. We parametrize the setup with a coupling of T = 0.3, 
a charging energy U = 3 (lOT) and a half-bandwidth 
D = 5 (15r) similar to Ref. [95l Band cutoff effects at 
D > U are expected not to play a role in the dynamics 

muz]- 


0.35 
0.3 
0.25 
„ 0.2 
0.15 
0.1 
0.05 
0 

t t 

FIG. 2. Current as a function of time t after a voltage quench, 
for temperatures T indicated. Parameters U = 3, = 0, 

r = 0.3, D — 5 and V = 2. (a) Monte Carlo (solid symbols, 
full lines), NCA (dotted lines) and OCA (dash dotted lines) 
results, (b) Zoom at times t = 0.3 — 0.9. (c) Behavior of the 
NCA and OCA currents up to times t = 5. 

Fig.§ shows the behavior of the current / after a volt¬ 
age quench. In panel (a) we show / as a function of time t 
following a voltage quench from R = 0 to R = 2 for a set 
of temperatures 0.07 < T < 1 (between 0.7 and 11 K). 
Results from semi-analytical approximations are shown 
as fine dotted lines (NCA) and dash-dotted lines (OCA) 



for comparison [Ml El- We observe that for all temper¬ 
atures the current equilibrates at time t « 0.6 (« 2.5 ps) 
well within reachable times of ~ 1 (4 ps). As temper¬ 
ature is lowered from T = 1 (11 K) to T ^ 0.2 (2 K) 
the current at fixed times and in steady state increases. 
Further reducing T by a factor of five yields no addi¬ 
tional increase of current, illustrated by panel (b). We 
attribute the fast equilibration time to the high Kondo 
temperature of this mixed valence system m 

Results from NCA and OCA correctly capture the 
short-time behavior and the overall shape of the cur¬ 
rent but underestimate both the transient and the steady 
state value by « 20%. This is known for systems with 
U < D [M], although the quality of these methods 
will presumably improve in the strong interaction limit 
and by including vertex corrections in OCA [93] . Both 
results are shown in panel (c) for times much longer 
than presently accessible by QMC. No additional time- 
dependence is visible, illustrating that our calculations 
are able to reach steady state. 



FIG. 3. (a) Transient current after a voltage quench at times 
t = 0.2, 0.4,0.6,0.8 and t = 1.0 as a function of temperature 
at 1/ = 3, = 0, F = 0.3, D = 5 and voltage 0 = 2. (b) 

Normalized current I{T)/I{T = 0) and equilibrium magnetic 
susceptibility x(r)/x(T = 0) as a function of temperature. 
Dashed line shows the value of 1/2 and vertical arrows show 
the crossing points of this line. 


The temperature dependence of the current is ana¬ 
lyzed in Fig. [^ The saturation of the low-T current is 
clearly visible in panel (a), where different traces show 
the transient current obtained at different times. We 
identify the saturation of the current in the low-T regime 
T = 0.07 — 0.2 with Kondo behavior and the temperature 
at which the current reaches half of its saturated low-T 
value as the Kondo temperature T/y [4|[I6|. 

Fig. i) shows the temperature dependence of the 
current ///(O) and equilibrium magnetic susceptibility 
x/x(0) normalized to the respective zero-temperature 
values. The magnetic susceptibility, defined as a re¬ 
sponse to the infinitesimal local magnetic field h, x = 
saturates at low temperatures [ini[ii] , sim- 

h —^0 

ilarly to the current. For K = 2 all normalized current 
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values for different times collapse on a single curve. This 
shows that Kondo behavior can be detected based purely 
on short-time transient dynamics. The Kondo tempera¬ 
ture Tr- determined by the current is time-independent 
and occurs at around twice the value of as defined 
from the point where the equilibrium magnetic suscep¬ 
tibility reaches half of its zero-T value, as expected in a 
mixed-valence regime [IlllMlISz!. 
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FIG. 4. (a) Quantum dot charge and (b) current aX U = 

3,ed = 0.0, r = 0.3, D = 5.0 and temperature T — 0.1 for 
voltages V = 0.8, 2,4, 6 (equal to 2U), 8, and 10 (equal to 
2D). 


We proceed with studying the impact of the magni¬ 
tude of the applied bias voltage. Fig. shows the time 
dependence of the quantum dot occupancy (a) and cur¬ 
rent (b) at a range of voltages between 0.8 and 10. At 
voltages V < 2 (linear response) the occupancy retains 
its equilibrium value and the current shows a monotonic 
rise and saturation to the steady state value, which in¬ 
creases with the applied voltage. Larger voltages V > 2 
demonstrate linear response behavior only at small times 
t < 0.3. At larger times the nonlinear behavior is visible: 
the current decreases to the steady state value, which be¬ 
comes smaller with applied F, illustrating the breakdown 
of conductance at large voltages [M] . 

The temperature dependence of the current outside of 
the linear response regime is analyzed in Fig.j^ Panel (a) 
shows I/I{0) at F = 6 as a function of temperature T. 
Different traces show different transient times, between 
0.2 (1 ps) and 1 (4 ps). Similarly to the linear response 
case (Fig. |^, the current exhibits the low temperature 
saturation at all times. The Kondo effect is therefore vis¬ 
ible in transient dynamics for all voltages. Outside linear 
response, the temperature at which the current saturates 
is strongly time-dependent. We denote the character¬ 
istic saturation temperature as Tt and define as before 
I{Tt,t) = I[T = 0,t)/2. At short times (panel (b)) Tt 
is the Kondo temperature Tk for all voltages. As time 
increases Tt increases by « 60% and reaches its steady 
state value of Tt ~ 2.9 for F = 6 (c). Further increase 
of the voltage results in a non-monotonic temperature 
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FIG. 5. (a) Normalized current I/I{T = 0) as a function of 

temperature at times t = 0.2,0.4,0.6, 0.8,1.0 and U = 3, Ed = 
0.0, F = 0.3, D = 5.0 and voltage F = 6.0. The dashed line 
indicates the value of 1/2. (b) Normalized current I/I{T = 0) 
as a function of temperature at t = 0.2 and (c) t = 0.8 for 
different voltages. 


dependence of the current |39j (not shown). 
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FIG. 6. (a) Current saturation temperature Tt {T at half 

of zero-temperature current value) as a function of time at 
different voltages at U = 3, = 0, F = 0.3, D = 5. Solid 

lines: estimates of the steady-state Tt. (b) Tt as a function of 
voltage at a set of times. Solid line: linear response Tt = Tk, 
dashed line: linear fit at F > 2. 

In Fig.j^we show the time (a) and voltage (b) evolution 
of Tt. For small voltages, there is no time-dependence of 
Tt. In contrast, as F is increased up to 2U, Tt rapidly 
rises and equilibrates. Estimates of the steady-state Tt 
as defined by the value at which the current reaches half 
of its low-T value for F = 2,4,6 are given by the hor¬ 
izontal lines. Fig. shows a separation of short- and 
long-time behavior: no change in Tt is observed for short 
times t < 0.4 (1.5 ps), whereas Tt in the steady state is 
found to be increasing with F for F > 2 in agreement 
with experiments |38j and predictions from the real time 
renormalization group |39| . 

Conclusions. We have described the transient dynam¬ 
ics of a quantum dot in the mixed valence regime follow¬ 
ing the instantaneous application of a bias voltage in a 
range of temperatures below and above the Kondo tem¬ 
perature Tk ■ We have observed the full dynamics of the 
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system from equilibrium to steady state. At all times 
and all voltages below the lead bandwidth, the current 
saturates at low temperature, exhibiting Kondo behav¬ 
ior. In linear response the saturation temperature Tj at 
which the current reaches half of its zero-temperature 
value is the Kondo temperature Tk- Outside of the lin¬ 
ear response regime T* has a strong time dependence and 
connects the equilibrium Kondo temperature to the in¬ 
creased steady state value [551ISH] ■ The temperature T( 
describes the interplay of non-equilibrium (applied volt¬ 
age) and strong correlations in the system. 

The results presented here are exact within the error- 
bars and the current and time-dependent Tt are directly 
accessible in time-resolved experiments (HI US] . The dy¬ 
namics of the mixed valence system is fast - the equili¬ 
bration occurs on a picosecond scale. We believe that the 
same physics can be observed at much larger time scales 
by lowering the level spacing ea and, consequently, en¬ 
tering the Kondo regime, thereby decreasing Tk and ex¬ 
ponentially increasing the relaxation time r oc exp(— Tr-) 
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